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Abstract. We discuss an NNLO realization of the general mass scheme S-ACOT-^ for treatment 
of heavy-flavour production in neutral current deep-inelastic scattering. Practical implementation of 
the NNLO calculation is illustrated on the example of semi-inclusive structure functions F2 C (x,Q) 
and Fic(x,Q). 
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Introduction. Correct computation of heavy-quark contributions to deep-inelastic 
scattering in the global PDF analysis is essential for predicting precision cross sections 
for W and Z boson production at the LHC [1]. 

Quark mass effects on DIS cross sections are comparable to next-to-next-to leading 
order (NNLO) contributions, therefore they must be included consistently in perturba- 
tive computations. The S-ACOT-^ scheme [2, 3, 4] is the default general-mass (GM) 
framework of CTEQ global PDF analyses; it is an improved formulation of the original 
ACOT scheme introduced in [5, 6]. 

S-ACOT-^ is an algorithmic procedure to account for heavy-quark masses in DIS 
structure functions that descends directly from the proof of QCD factorization for DIS 
with massive quarks [2]. Compared to the earlier ACOT and S-ACOT schemes, the 
S-ACOT-^ scheme adds a phenomenologically important requirement that all heavy- 
quark scattering contributions (including those with collinear approximation for heavy 
quarks) satisfy energy-momentum conservation near the quark production threshold. 
In Ref. [7], we show that the constraints from energy-momentum conservation can be 
included directly in the QCD factorization theorem. Thus, the S-ACOT-^ scheme is 
proved to all orders on the same footing as the original ACOT scheme. 

We also present, for the first time, S-ACOT-^ numerical predictions of NNLO accu- 
racy for heavy-quark neutral-current DIS structure functions, F2 C (x, Q) and Fl c (x, Q 2 ). 

S-ACOT-^ is attractive for phenomenological applications because of its relative sim- 
plicity compared to other schemes like BMSN [8] or TR [9]. It achieves a great sim- 
plification by replacing coefficient functions with incoming heavy-quark lines by their 
zero-mass (ZM) expressions, while requiring that the light-cone momentum fraction £ 
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in the PDFs is always within the interval % < t, < 1 that satisfies energy- 

momentum conservation [4, 10]. Here and in the following, a rescaled "Bjorken-x" vari- 
able % = x (1 + 4m^/ Q 2 ) is introduced. This arrangement allows to include all the phys- 
ical scattering channels, while excluding kinematically prohibited £ values. 

The S-ACOT-^ scheme assumes one value of the number of active quark flavors, Nf, 
and one PDF set in each Q bin, in contrast to other possible GM schemes. Matching of 
massive 4-flavor predictions on massive 3-flavor predictions at low Q, and on massless 
4-flavor predictions at large Q is realized in S-ACOT-^ through cancellations among 
certain classes of Feynman diagrams. These cancellations are improved by using the 
optimal rescaling variable %, but also hold for other forms of the rescaling variable. 
Practical implementation is easy to handle, and all NNLO Wilson coefficient functions 
that are needed are available from literature. 

Summary of the computation. If the components of inclusive F2±{x, Q 2 ) = F are 

classified according to quark couplings to the photon [1 1], F can be written as a sum of 
light and heavy-quark coupled contributions, 

N, 

F = ^F l +F h , (1) 
l=l 

where 

Fl = e 2 i £ [Q, a ® fa/p] (x, Q), F h = e\ £ [C M <g> f a/p ] (x, Q) . (2) 

a a 

Here Q ^ are Wilson coefficient functions, and f a / p are PDFs. Perturbative expansion up 
to G(a 2 ) leads to 

Fl 2) = 4 {% h {2) ® (fh/p + fh/ p ) +Cf) ® E + cg ®f g/p ] , (3) 

where lowercase c^ 2 ] and uppercase C^ a - h represent ZM coefficient functions and mas- 
sive coefficient functions, constructed from results inRefs. [12, 13, 14] and Refs. [8, 15, 
16, 17], respectively. E is the singlet PDF combination. A detailed derivation of these 
coefficient functions is presented in Ref. [7]. A similar expression for F\ is given by 

F« 2 » = 4 {f,,, +fl/p ) + 8 1 + c g> S /„, } . (4) 

where cf^ g , c PS '^ and cf^^ are constructed from components available in Refs. [12, 
13,14,7]. 

Numerical results. In the following, we show representative NNLO predictions for 
F2c and Fl c , computed using Les Houches toy PDFs [18, 19] that are evolved with 4 
active flavors by HOPPET computer code [20]. In Fig.l, NNLO S-ACOT-^ predictions 
for Fic (left panel) and Fu (right panel) are shown vs. Q by blue solid lines. ZM 4-flavor 
predictions at NNLO are shown by purple long-dashed lines, and 3-flavor (fixed- flavor- 
number, or FFN) predictions by red short-dashed lines. Fig.l shows that S-ACOT-^ 




FIGURE 1. Upper insets: semi-inclusive F2 C (x,Q) and Fi c (x,Q) at NNLO as a function of Q at 
x = 10~ 2 . Lower insets: ratios of ZM and FFN predictions to the respective S-ACOT-;^ predictions. 



smoothly interpolates in between the FFN and ZM schemes, and that it reduces to the 
FFN scheme at Q « m c and to the ZM scheme when Q^> m c . Lower insets show ratios 
of the FFN and ZM predictions to the respective S-ACOT-^ predictions, to elucidate 
differences in the intermediate region. 

Fig. 2 shows a remarkable reduction of the factorization scale dependence of i^Lc- 
Without extra tuning of the factorization scale, the S-ACOT-^ prediction is close to 
FFN and other NNLO schemes at Q w m c . Here scattered symbols in the left panel 
correspond to predictions based on MSTW/TR' NNLO coefficient functions (sea-green 
triangles) [21] and FONLL-C NNLO coefficient functions (blue circles) [11]. The solid 
black line is the S-ACOT-^ NNLO prediction corresponding to a reference factorization 
scale /i = a/ Q 2 + m 2 . The green band is the theoretical uncertainty in the S-ACOT-^ 
prediction due to the variation of the scale in the range Q < ji < \J Q 2 + 4m 2 . The purple 
band represents scale variations in the FFN prediction at NNLO around the reference 
scale values indicated by the magenta short-dashed line. The light blue band around the 
blue-dashed line represents the S-ACOT-^ prediction at NLO and its scale uncertainty. 

Conclusions. We have shown that the S-ACOT-^ scheme is fully compatible with 
the QCD factorization theorem, and that an NNLO computation of Fi^Lc in the S-ACOT- 
X scheme is viable. S-ACOT-^ predictions at NNLO are stable and show significant 
reduction in the factorization scale dependence (see Fig. 2), compared to NLO computa- 
tions. This is the most challenging component of the CTEQ global analyses at NNLO. 
A full description of NNLO S-ACOT-^ computations is available in [7]. 
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FIGURE 2. Factorization scale dependence for semi-inclusive F2,lc(x,Q) at NNLO as a function of x 
at Q = 2 GeV. 
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